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We report a measurement of the production cross section for b hadrons in pp collisions at = 
1.96 TeV. Using a data sample derived from an integrated luminosity of 83pb~^ collected with the 
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upgraded Collider Detector (CDF 11) at the Fermilab Tevatron, we analyze b hadrons, Ht, partially 
reconstructed in the semileptonic decay mode Hb ^~ D°X. Our measurement of the inclusive 
production cross section for h hadrons with transverse momentum pr > 9 GeV / c and rapidity 
\y\ < 0.6 is o" = 1.30 /ib ± 0.05 /^b(stat) ± 0.14/ib(syst) ± 0.07/ib(C), where the uncertainties are 
statistical, systematic, and from branching fractions, respectively. The differential cross sections 
dcr/dpT are found to be in good agreement with recent measurements of the Hb cross section and 
well described by fixed-order next-to-leading logarithm predictions. 

PACS numbers: 14.40.Nd, 12.15.Ff, 12.15.Hh, 13.20.He 

I. INTRODUCTION 

Measurements of the 6-hadron production cross section in high energy pp collisions provide an excellent test of 
perturbative quantum chromodynamics (QCD). Next-to-leading order (NLO) calculations of the &-hadron cross section 
[l|, 0] have been available for more than a decade. These NLO calculations have been recently supplemented with 
corrections for terms proportional to a1[as\og{pT/^b)]^ and aj[as log(pT/mf,)]'^ with fc > 1, where is the strong 
coupling constant, pt is the transverse momentum of the b hadron, and mi, is the mass of the bottom quark. These 
calculations are referred to as fixed order next- to- leading logarithm (FONLL) 0]. 

A number of 5-hadron cross-section measurements has been performed at the Tevatron. In the 1992-1996 running 
period, Tevatron Run I, measurements were performed at ^/s — 1.8 TeV by the CDF and DO experiments using 
several different techniques and final states [1, [a, B 13 • These measurements were consistently higher than theoretical 
expectations from NLO calculations. 

Since 2001, the Tevatron has been running at = 1.96 TeV. To date, the CDF Collaboration has made two 
inclusive 5-hadron cross-section measurements at the higher center-of-mass energy. The first was an inclusive mea- 
surement using —> J/ipX, J/ip — > ^JL^ decays, where Hf, denotes a h hadron [8|. This measurement was performed 
for \y\ < 0.6 and was the first to map out the 6-hadron cross section down to zero transverse momentum. The 
Hb — > J/ipX cross-section result, along with an improved FONLL calculation fsj, prompted studies [lOjIiii suggesting 
that prior discrepancies had been resolved with improved measurement and calculational techniques. 

More recently, the CDF Collaboration has performed a measurement of the meson cross section using 740 pb~^ 
of data through the fully reconstructed decay chain i?+ — > J/tpK^, with J/ijj — > M"*"/^", for transverse momentum, 
Pt{B^), greater than 6GeV/c and |y| < 1. With large statistics and a high purity signal, the total uncertainty on 
this measurement is smaller than 10% [T5]. The result is also in good agreement with FONLL calculations. The 
difference between the original NLO calculation and the current FONLL calculation arises from a number of different 
factors, as discussed in Ref. ilO]. 

While improved theoretical calculations compare favorably to recent measurements, more experimental input is 
needed. In particular, cross-section measurements using different Hi, decay modes can be complementary to one 
another since they might be acquired using different trigger paths and be subject to different background contributions. 
For example, Ref. [13[ showed that previous 6-hadron cross-section measurements were inconsistent when comparing 
semileptonic and J/^p final states. 

In this paper, we present the first measurement of the Ht cross section using semileptonic decays at ^/s = 1.96 TeV. 
This analysis takes advantage of the distinct semileptonic 6-hadron decay signature and provides a measurement that 
is complementary to those made in H}, J/ipX modes. We reconstruct signals of semileptonic 6-hadron decays, 
Hb fj,~DX, where charm mesons, and D*'^ , are reconstructed in fully hadronic modes, using a data sample 
acquired with a dedicated semileptonic trigger path. Throughout this paper, any reference to a specific charge state 
implies the antiparticle state as well. 

Previous measurements of the 5-hadron cross section at the Tevatron using semileptonic decays were performed using 
inclusive Hb ix~ X 0,0] and the semielectronic Hb e~D^X [5] modes. There are no previously published 6-hadron 
cross-section results from the Tevatron utilizing the Hb — > D'^X mode. We improve upon prior semileptonic results 
with a comparable data sample by extending the measurement to lower values of pxiHb) than have been accessible 
in previous analyses. The trigger path used in this analysis, which will be described in detail below, accepts muons 
with px > 4GeV/c while previous measurements were limited to leptons with px > 8GeV/c. 

We begin in Sec. [Til with an overview of the measurement technique. The CDF II detector and trigger system will 
be described in Sec. IIIIl followed by a discussion of trigger selection in Sec. IIVI Section |V] will describe the event 
selection and candidate reconstruction. Sections IVII and IVIII then describe the measurements of signal reconstruction 
efficiency and acceptance, followed by an assessment of the background contributions in Sec. IVIIII In Sec. IIX| we 
present our results along with a comparison to theoretical predictions. 
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II. OVERVIEW OF THE MEASUREMENT 

We select events with an identified muon in conjunction with a fully reconstructed charm meson {D'^ or D*^ , 
generically referred to as a D meson) such that the ii~D topology is consistent with a 6-hadron decay. To measure 
the cross section, we count the number of signal events and correct for detector acceptance, identification inefhciency, 
and contributions from background sources. Incorporating these factors, along with the measured pp luminosity, C, 
we extract the cross section, cr, times branching ratio (B) in the following way: 

a{pp^mX)xB= ^^\-/^^'^ , (1) 

where N is the number of observed candidate n^D^ (fi^D*^) events and fbck is the background fraction, i.e. the 
fraction of events not originating from the decay of a & hadron. The acceptance a is defined as the fraction of events 
in which the muon and the charm decay particles pass through active regions of the CDF II detector. The detection 
efficiency e is the fraction of these events that are reconstructed. 

The product of branching ratios, B, is shorthand for B{Ht, — > ijl~D^X) x B{D'^ — > K~-k^) in the nT mode and 
B{Hi, -> ^ji'D*+X) X B{D*+ D^7r+) x B{D° K-n+) for the mode. The generic b hadron, Hb, is an 

admixture of all weakly decaying b hadrons, including i?", B^, B^, and Ab. To extract the Hb cross section, we 
utilize the measurement of the inclusive Hb fji" D'^X branching ratio from the DELPHI and OPAL experiments 
[31 . To apply this branching ratio to Tevatron data, we assume that the fragmentation fractions at high energy 
{pT{Hb) > 9 GeV/c) are the same at LEP and the Tevatron. Charm branching ratios are taken from CLEG, ARGUS, 
and ALEPH measurements [H, [l^ . 

Under the assumption that b and b quarks are produced in equal rates at the Tevatron, we include both /i+Z) and 
fi~D° final states and introduce the factor of 2 in the denominator of Eq. ([T]) to report the cross section for a single 
6-hadron flavor. 

To extract the number of candidate signal events, N, we reconstruct Hb fi^D'^X with D'^ if^7r+, as well 
as Hb ji^ D*^ X with D*^ —> D'^tt^ and D'^ K^tt^ , events using data that were acquired via a dedicated 
semimuonic Hb trigger that requires a well measured muon and a track displaced from the primary vertex. The 
D*^ sample is a subset of the K~tt~^ sample. We use these two samples to make separate, but correlated, 

measurements . 

The selection criteria are chosen to preferentially select Hb decays and reject combinatoric backgrounds. There are 
two types of physics backgrounds that exhibit a. ^~ D signature that is similar to that of the expected signal. Gne 
arises from direct pp — > ccX production where one charm quark decays to a muon and the other fragments into a 
D meson. The second physics background arises from pp — > bbX events where one b hadron decays to a muon while 
the other decays to a D meson. The contributions from these sources are assessed using a combination of data and 
Monte Carlo (MC) simulation. 

The acceptance a is determined by using MC simulation. We generate b quarks based on input distributions 
in transverse momentum pT and rapidity y taken from theoretical calculations as well as previous cross-section 
measurements. Fragmentation and decay are performed as part of the MC simulation. The passage of generated 
decay products in the CDF II detector is then simulated. Based upon the output of these simulated events, the 
fraction of produced events that fall within the active detector volume is calculated. 

We factorize the determination of the trigger and detector efficiency, e, into a product of relative efficiencies and an 
absolute efficiency. The advantage of this technique is that we can measure relative efficiencies from the data (as will 
be described in Sect. IVip and minimize our reliance on simulation. 

The sensitivity of this analysis is limited by the systematic uncertainties. To limit our dependence on theoretical 
input and on detector simulation we estimate systematic uncertainties from data wherever possible. 

III. THE CDF II DETECTOR 

The CDF II detector comprises a charged-particle tracker in the magnetic field of a solenoid surrounded by a full 
coverage calorimeter and by muon detectors. In this section we give a short description of those detector components 
that are relevant to this analysis. A more detailed description of them and of the entire detector is given in Refs. (TtI . 
[18, 19, 20, 21]. 

We utilize a cylindrical coordinate system with the z axis oriented in the direction of the proton beam and the r-t/f 
plane perpendicular to the beamline. We also define a polar angle, 9, measured from the z axis and pseudorapidity 
is defined as 77 = — ln[tan(0/2)]. 
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Two devices inside the 1.4 T solenoidal magnetic field are used for measuring the momentum of charged particles: 
the silicon vertex detector (SVX II) and the central tracking chamber (COT). The SVX II consists of double-sided 
microstrip sensors arranged in five cylindrical shells with radii between 2.5 and 10.6 cm. The SVX II detector is divided 
into three contiguous sections along the beam direction for a total z coverage of 90 cm. The COT is a cylindrical 
drift chamber containing 96 sense wire layers grouped into eight alternating superlayers of axial and stereo wires. The 
active volume covers \z\ < 155 cm and 43 to 132cm in radius. 

The central muon detector (CMU) surrounds the central electromagnetic and hadronic calorimeters. The calorime- 
ters have a depth corresponding to 5.5 interaction lengths. The CMU detector covers a pseudorapidity range |r/| < 0.6 
and is segmented into two barrels of 24 modules, each module covering 12.6° in cj). Every module is further segmented 
into three submodules, each covering 4.2° in cj) and consisting of four layers of drift chambers. The smallest drift unit, 
called a stack, covers 1.2° in (p. Adjacent pairs of stacks are combined together into a tower. A track segment is 
identified by at least two hits out of the four layers of a stack. A second set of muon drift chambers, the CMP, covers 
the same 77-range as the CMU and is located behind an additional steel absorber of ~3.3 interaction lengths. Muons 
that produce a track segment in both the CMU and CMP systems are called CMUP muons. A third set of muon 
drift chambers, the CMX, covers a pseudorapidity range 0.6 < \ri\ < 1.0 and is separated by 6.2 interaction lengths 
from the nominal interaction point. 

Luminosity is measured using gas Cherenkov counters (CLC) mounted at small angles to the beamline. The 
CLC measures the rate of inelastic pp collisions. The inelastic pp cross section at ^/s — 1960 GeV is scaled from 
measurements at ^/s = 1800 GeV using the calculations in Ref. [22|. The integrated luminosity is determined with 
negligible statistical uncertainty and a 6% systematic accuracy [23] . 

The CDF experiment uses a three- level trigger system. At level 1 (LI), data from every beam crossing are stored 
in a pipeline capable of buffering data from 14 beam crossings. The level 1 trigger either rejects events or copies them 
onto one of the four level 2 (L2) buffers. Events that pass the level 1 and level 2 selection criteria are sent to the 
level 3 (L3) trigger, a cluster of computers running a speed-optimized version of the offline event reconstruction code. 

To select heavy flavor events, we rely heavily upon charged-particle tracking in the trigger. At level 1, charged 
tracks with pr > 1.5GeV/c and \ri\ < 1 are found by the extremely fast tracker (XFT) 24], which uses information 
from the axial wires in the COT to perform r-0 track finding with high efficiency (> 90%) and good transverse 
momentum resolution [Spx/Px — 0.017). The track extrapolation system (XTRP) takes the track information from 
the XFT and provides extrapolation information so that the XFT tracks can be matched to track segments found 
in the muon detectors 25]. Tracks are matched to track segments in the CMU, CMP and CMX to identify muon 
candidates at level 1. 

At level 2, the Silicon Vertex Trigger (SVT) combines information from the axial layers of the SVX II with the 
charged-particle tracking information from the XFT. Owing to the high precision information from the silicon detector, 
the SVT provides a precise measurement of the track impact parameter do (distance of closest approach to the 
beamline) in the transverse plane with a resolution of approximately 50 /xm. This resolution arises from approximately 
equal contributions of 35 /.tm from resolution and intrinsic spread of the position of the pp interaction. Tracks with 
large impact parameter are utilized to identify heavy flavor decay and strongly reject light (u, d, s) quark events. 

For the data sample utilized in this analysis, level 3 event reconstruction included full COT tracking but did not 
include tracking in the silicon detectors. Events accepted by the level 3 trigger are stored for subsequent analysis. 

The flexibility of the CDF II trigger permits the selection of samples with no leptons (all hadronic modes), one 
lepton (semileptonic modes) and two leptons {J/ip — *■ ^J■^ ^J'~ and T —^ n'^ ii~ events.) For this analysis we utilize all 
three types of triggers, and we describe them in the following section. 

IV. TRIGGER PATHS 

A single trigger path is defined as a specific set of level 1, level 2, and level 3 selection criteria in the CDF II trigger 
system. The primary trigger path utilized in this analysis is referred to as the /iSVT path. Data acquired through 
other trigger paths are utilized in this analysis to determine the selection efficiencies described in Sect. IVII 

/iSVT Path. At level 1, this path identifies a muon candidate with px > 4GeV/c by requiring that an XFT track 
be matched to segments in the CMU and CMP detectors. At level 2 the trigger utilizes information from the SVT 
to identify displaced tracks. We require a single track (which is not the muon candidate) with px > "2 GeV/c and an 
impact parameter between 120 fim and 1 mm to be identified, where the impact parameter measured by the SVT is 
relative to the beamline as determined over a large sampling of events. At level 3, full event reconstruction is done, 
except for track reconstruction in SVX II. In order for the event to be accepted, level 3 requires confirmation of the 
muon identified at level 1 and confirmation of the XFT track identified at level 2. The invariant mass of the two 
tracks must be less than 5GeV/c^, the azimuthal opening angle between the two tracks, Acj), must be 2° < Acf) < 
90°, and the tracks must be consistent with coming from the same pp interaction vertex. 
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Inclusive Muon Path. This trigger requires that one muon candidate with pT > 8GeV/c be identified by the 
XFT, CMU, and CMP. This path provides a sample oi J/tp ^ decays where only one of the muons is required 

in the trigger. The second muon in the event is unbiased, and can be utilized to directly measure the efficiency of the 
trigger selection. 

Heavy Flavor Path. This trigger path is for hadronic decays of Hi, and D hadrons. Two oppositely charged XFT 
tracks with pr > 2 GcV/ c and Aip < 135° are required at level 1. Those tracks are then required to have 120 pm < 
do < 1 mm as measured by the SVT at level 2. The level 3 selection requires confirmation of the level 2 quantities 
using full reconstruction. From this path, we select large samples of D'^ K'tt^tt'^, (/)7r+ K~^K~7r'^ , and 

Df — > 07r+ X+if ^7r+ decays. Since only two of the three tracks are required in the trigger, the third track is 
unbiased with respect to the trigger. We can therefore use this sample to measure the efficiency of the SVT. 

Hi, Semileptonic Backup Path. This trigger path is used to determine the level 3 selection efficiency. Its 
requirements are the same as the /iSVT trigger, except that it has no requirements at level 3, and only a fraction of 
events that pass the level 1 and level 2 trigger requirements are recorded. 

J/tp Dimuon Path. The J/tp trigger paths which require one muon in the CMU system and a second muon 
in either the CMU or CMX system are used to measure the SVX II offline efficiency. The J/ip CMU-CMU trigger 
requires, at level 1, at least two CMU muons be found with px > 1.5GeV/c. There are no additional selection 
criteria applied at level 2. At level 3, an opposite sign muon pair is required with invariant mass between 2.7GeV/c^ 
and 4.0GeV/c^, which is efficient for J/ip and ip{2s) decays into dimuons. The two muon tracks must also have 
\Azo\ < 5cm and A0o < 130°, where Azq is the separation of the two tracks along the bcamlinc (r — 0) and AcpQ is 
the separation of the two tracks in azimuth at r = 0. The J/ip CMU-CMX trigger is similar, except that instead of 
two CMU muons, one CMU and one CMX muon must be identified. The CMX muon must have pt > 2.0GeV/c. 

V. EVENT SELECTION 

This analysis uses data derived from 83 pb"-'^ of pp collisions at ^/s = 1.96 TeV that were collected between October 
2002 and May 2003. Thanks to the trigger utilized in this analysis, this relatively small sample of data is sufficient 
to significantly improve upon previous cross-section measurements in semileptonic channels, and yields a result that 
is limited by the systematic uncertainties. 

We begin the signal reconstruction by requiring the candidate muon, pion, and kaon tracks to satisfy COT fiducial 
and quality selection criteria, which include requirements on the number of hits assigned to the track and a loose 
selection on fit quality of the track. Muon selection criteria are motivated by the trigger requirements. The muon 
track must have pt > 4GeV/c and |?7| < 0.6 and must be matched to CMU and CMP track segments that satisfy 
fiducial and quality selection criteria, which include requirements on the number of hits used in the track segment 
and the quality of the match between the central track and the CMP track segment [1^ . The candidate muon must 
also be matched to a CMUP muon found by the level 1 trigger. 

In order to reduce combinatoric background, the pion and the kaon candidates must have px > l.OGeV/c and 
\r]\ < 1.0, must satisfy SVX II fiducial and quality cuts, which include requirements on the number of hits assigned 
to the track and a loose selection on the fit quality of the track [2^ . The two tracks must both come from a common 
displaced vertex and be consistent with the decay of a charmed meson. In addition, at least one of the candidate tracks 
must be matched to an SVT track that passes the displaced-track-trigger requirements, have 120 /im < do < 850 /im, 
and have hits in all of the SVX II layers used by the SVT. The pion and kaon candidates as well as the pion and 
muon candidates must have opposite charges. We combine the pion and kaon candidates to form — > K~Tr~^ 
candidates. We form the combined momentum vector of the muon and the candidate and require that this vector 
have |7;(/Lt~D°)| < 0.6. The number of D° signal events is determined by fitting the X~7r+ invariant mass distribution 
between 1.74GeV/c^ and 1.98GeV/c^ to a Gaussian peak plus a hnear background. 

The mass plots are shown for different regions of pt{iJ'~ D^) in Fig. [1] While the signal yield is rather low in 
the highest pT region, clean signals are observed in the lower pT regions. The signal yields are reported in Table [H 
Figure [H shows the mass distribution for all events with pt{h'D^) > 9.0 GeV/c and < 0.6. 

For the decay D*^ D'^tt^, K^n^, the muon and reconstruction is the same as described above. We 

additionally require the presence of the pion from the D*~^ iD°7r+ decay with px > 400 MeV/c and \ri\ < 1.0 that 
passes COT fiducial and quality selection criteria. Since this pion tends to have low pt, we refer to this as the "soft 
pion" [ti so ft-) We require that the muon and soft pion have opposite charges and only consider events where the 
K~TT^ invariant mass is between 1.82 GeV/c^ and 1.90 GeV/c^. As can be seen in Fig. [1] this mass window is wide 
enough so that no real decays are lost. We require \Azq\ between any two of the decay products of the Hb must be 
less than 5 cm. In order to minimize the background, we look at the mass difference Am between the K^tt^tt'^^j^ and 

the K~7r'^. We fit the mass difference distribution to a Gaussian plus the background shape a^/Ajri^^rn^e^''^'^~™'^\ 
where a and b are free parameters in the fit to the data. The mass difference plots are shown in Fig. [3l In order 
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to make comparisons between the two measurements, we present the data in px bins of the system. The 

yields are included in Table [H Figure 2] shows the mass difference for all events with p7-(/i"Z?") > 9.0 GeV/c and 
\tj{h-D")\ < 0.6. 

VI. EFFICIENCY MEASUREMENTS 

The efficiencies of the CDF II trigger and detector components for the cross-section measurement are determined 
whenever possible by using data collected by complementary trigger paths. The trigger paths utilized are described 
in Section ITVl 

We have separated the overall efficiency, e, into nine separate measurements, denoted ei-eg. The first eight terms 
are extracted using relative efficiencies and the ninth term is measured as an absolute efficiency. To calculate es, we 
measure the efficiency for events to pass selection criterion eight relative to events that already have passed criterion 
nine. This yields es — Cs'^'eg- Repeating this procedure for selection criterion seven, we calculate the relative efficiency 
for events to pass selection seven relative to events that have already passed selections eight and nine: ej = Cy'^'eg. 
Repeating this process for each of the selection criterion, we arrive at a final efficiency which is 

e = e,f[erK (2) 

1 

The absolute efficiency measurement (eg) is the efficiency for finding a charged-particle track in the COT. This absolute 
efficiency is determined using a combination of data and MC simulation based upon GEANT [13] ■ A brief description of 
these efficiency measurements is provided below. Details of all of the efficiency measurements and parameterizations 
may be found in Ref. fz^ . We perform these measurements in bins of kinematic variables. For example, the level 1 
trigger efficiency depends upon the momentum of the muon. The efficiency corrections are then applied according to 
the kinematics of each event. 

A. Level 1 Efficiency (ei) 

For the /xSVT trigger path, the level 1 muon trigger requirement consists of an XFT track with pr > 4GeV/c that 
is matched to track segments in the CMU and CMP subdetectors. We measure the efficiency of this trigger with 
respect to the offline CMUP muon reconstruction efficiency. To do this, we reconstruct the J/ip M^A*^ signal in 
events that were collected via the inclusive muon trigger path. The trigger muon [px > 8GeV/c) in these events is 
biased for the efficiency measurement, but the other decay muon of the J/V' can be used to determine the trigger 
efficiency. We refer to the J/tp decay muon that is unbiased by the trigger as the "probe" muon, since it can be used 
to directly measure the trigger efficiency. We then compare the number of times the probe muon satisfied the trigger 
to the number of probe muons that were within the acceptance of the trigger. 

We require both muons of the J/ip to pass COT, CMU, and CMP fiducial and quality selection criteria. The 
probe muon must also have \ri\ < 0.6 and pr > ^ GeV/c. We find a signal yield of approximately 1900 events. 
Sideband subtraction is used to parameterize the efficiency in terms of muon pT and rj. We find the plateau efficiency 
{PT > 5GeV/c) to be 90%. 

B. Level 2 XFT Efficiency (£2) 

For the /iSVT signal path, the XFT tracks found in level 1 are used as input for the level 2 SVT trigger. The 
trigger path requires at least one XFT track, not associated with the muon candidate, having > 2 GeV/c; the SVT 
further requires that the impact parameter for this track be between 120// m and 1 mm. The overall efficiency of this 
selection is given by the product of the XFT and SVT efficiency. We require that either the kaon or pion from the 
-D" decay satisfy this selection. 

The XFT efficiency is measured for pious and kaons with pT > 2.0 GeV/c and |ry| < 1.0. Given their different 
ionization energy losses, we treat kaons and pions separately in this measurement. In order to obtain pure samples 
of kaons and pions, we reconstruct — > K~tt~^tt^ events, D'^ 'P'^'^i 4' ~^ K^K~ events, and Z?+ — > (/>7r+, 
(j) K~^K~ events that were collected by the heavy flavor trigger path. We have roughly 200 000, 9500, and 18 
000 signal events respectively in the three modes. Two of the tracks from each decay are required to pass the trigger 
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requirements, and the third track is used to probe the XFT efficiency. We apply both COT and SVX II fiducial and 
quality selection criteria on the tracks that fired the trigger, and COT fiducial and quality selection criteria on the 
tracks used to probe the XFT efficiency. For the resonant decays, a cut is placed at ztlOMeV around the (j) mass. 
Sideband subtraction is used to measure the pion and kaon yields. We parameterize the efficiencies in terms of pr, 
•q and (j). We additionally parameterize the XFT cfhciencies over the time that the data was acquired. The average 
efficiencies are approximately 90% for pions and 85% for kaons. 

C. Level 2 SVT Efficiency (ea) 

The second component of the level 2 trigger efficiency is the SVT efficiency for tracks that have already been 
identified by the XFT. The SVT efficiency does not depend on the particle species, so we measure the efficiency using 
a high statistics sample of J/V' — > A*^A*~ events that were collected by the inclusive pt > 8 GeV/c muon trigger path. 
The technique here is similar to the technique utilized to measure the level 1 efficiency. 

We use tracks that have pt > 2.0 GeV/c, \rj\ < 1.0, and 120 fim < \do\ < 850 /im. The muon from the J/4> that 
did not satisfy the trigger is used as the probe track. On events where both muons satisfy the pr > 8 GeV/c inclusive 
muon trigger, both muons are used as probe tracks. We impose COT and SVX II fiducial and quality selection criteria 
on the probe track, and require the track to be matched to an XFT track. Therefore the SVT efficiency is measured 
relative to a track found offline with the required number of SVX II hits. There are approximately 71 000 tracks that 
pass these selection criteria. The efficiency is parameterized in terms of pt, ^Oj a-nd track isolation, which is a measure 
of the number of tracks found within an angular region of Scj) < 5°. For tracks with pT > 3GeV/c, the plateau 
efficiency for the SVT is approximately 90%. The rj dependence of the efhciency is small, and further mitigated 
because the tracks used in the efficiency measurement sample the detector in a manner similar to the signal. Since 
the SVX II hit efficiency is not part of the level 2 SVT efficiency calculation (it contributes to the SVX II efficiency 
described below) the inefficiency here is dominated by the track-finding and fitting algorithms utilized in the SVT. 

The track isolation dependence of the SVT efficiency requires special treatment. The occupancy in the silicon 
detector arises from a number of sources, including very low momentum tracks from the underlying pp interaction 
along with a spectrum of tracks originating from multiple pp interactions. These effects are not modeled reliably in 
the Monte Carlo simulation, and so we reweight the data based on the isolation of the signal events. In addition to 
the raw yields, Table U also lists the yields after correction for the SVT efficiency. 

D. Level 3 Efficiency (£4) 

At level 3, our trigger path requires that a muon with pT > 4GeV/c be reconstructed. In addition, a track identified 
in the COT (other than the muon) must be matched to an SVT track meeting the level 2 trigger requirements. The 
muon track and COT-SVT track must have an invariant mass of less than 5.0GeV/c^ and 2.0° < A0o < 90°. To 
find the efficiency of this trigger, we look at ^ 30 000 events that were acquired on the Hf, semileptonic backup 
path. We parameterize the efficiency in terms of the px of the muon and the px of the SVT track. For muons with 
Pt > 4.5GeV/c, the level 3 efficiency is approximately 97%. 

E. SVX II Efficiency for Triggered Pions and Kaons (es) 

In order to determine the efficiency for attaching SVX II hits to a COT track, we look at muon tracks from 
J/ijj fJ.~^ fi~ decays collected on the J/iJj dimuon trigger path. We use this sample in order to have a large sample of 
tracks with negligible fake rate. The tracks are required to pass COT quality and fiducial, as well as SVX II, selection 
criteria. We consider only tracks that are expected to pass through the four layers of the SVX II used in the SVT 
trigger. The average efficiency is observed to be 80%, where the inefficiency is dominated by dead silicon ladders and 
poor efficiency near the ends of the silicon ladders. To account for the dependence on SVX II coverage, we correct 
the efficiency as a function of the z position of the track at its origin. 

F. SVX II Efficiency for Untriggered Pions and Kaons (eg) 

We require hits in at least three layers of the SVX II on the decay track that is not required by the trigger. We 
find the efficiency by looking in our signal sample with and without the SVX II hit requirement for the second track. 
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This yields an SVX II efRciency for the untriggered track of (93.2 ± 1.0)%. The efficiency for the untriggered track is 
higher than that for the triggered track because of the fooser requirement on the number of SVX II hits. 



G. CMU efficiency (er) 

In order to find the efficiency for linking a muon COT track to a CMU track segment, we reconstruct J /ip M"*"/^" 
events using both the inclusive muon and semileptonic backup trigger paths. The muon tracks must pass COT fiducial 
and quality selection criteria. One of the decay muons from the J/iJj must satisfy the muon trigger requirements. For 
events acquired on the semileptonic backup trigger path, the probe track must satisfy the SVT trigger requirements. 
The probe tracks are extrapolated to the CMU, and are required to pass CMU fiducial selection criteria. The probe 
tracks must have pr > 4.0 GeV/c and \ri\ < 0.6. The efficiency is found to be independent of pr, 0, and rj, and so we 
use a single efficiency value of (79.5 ± 1.3)%. 



H. CMP efficiency (eg) 

We find the efficiency for linking a CMP track segment to a COT track that has already been linked to a CMU track 
segment by reconstructing J/ijj jJL^ n~ candidates in events collected by the J /ip dimuon trigger path. This trigger 
path does not require hits in the CMP, so either decay muon from the J /ip can be used to probe the CMP efficiency. 
We require the muon to pass COT and CMU fiducial and quality selection criteria. The muon is extrapolated out to 
the CMP and required to pass CMP fiducial selection criteria. We examine only muons with px > 4.0 GeV/c and 
|?7| < 0.6. The efficiency is parameterized in terms oi pt and ?/. The average CMP efficiency is 85%. 



I. COT efficiency (eg) 

To determine the absolute COT track-finding efficiency, we use a combination of data along with a detailed MC 
simulation of the CDF II tracking system. The techniques described here have been utilized in previous CDF analyses 
2^ . We have verified that the MC simulation accurate models the charged-particle kinematics of our signal sample 
26|. We have additionally performed detailed studies of the COT simulation and find that the data is well- modeled 
by the simulation [29l |. 

First, we use a MC simulation only technique to compare charged particles generated by the pythia MC program to 
charged tracks reconstructed after full MC detector simulation and event reconstruction. This yields our measurement 
of absolute COT track-finding efficiency. 

To assess the systematic uncertainty on this measurement, we use a combination of data and MC simulation to 
"embed" simulated MC kaons, pions and muons into data events acquired via the inclusive muon trigger path. The 
MC simulated COT hits from the simulated track are incorporated into the detector data from a real event. To 
accurately emulate track reconstruction in D'^ events, the embedded MC simulated track is placed within a region 
in rj-(t> space that is within -y^ Sjf + S(f>'^ < 1.4 of the trigger muon. The hybrid event, which consists of data plus 
hits from a single simulated track, is then reconstructed in the same way as our signal sample. The simulated track 
is considered to be correctly found if a reconstructed track is matched closely in px and (p to the simulated track 
embedded in the event. 

We vary the matching criteria to assess the rate at which tracks may be found but reconstructed with incorrect 
parameters. For tracks with px > l.OGeV/c (which includes all muons and kaons, and the pions from the decay of 
the D'^), the track-finding efficiencies derived from MC-only simulation and track embedding agree to within 1% and 
are insensitive to the details of the matching criteria. The COT efficiency for muons is approximately 99%. For pions 
and kaons with pr > 2 GeV/c the COT efficiency is greater than 90%. For pions with pr < 1.0 GeV/c, which are used 
in reconstructing the D*^ — > D^n^ decay, we observe a significant difference in efficiency estimates from the different 
techniques, and we use these variations to assess a systematic uncertainty on the low px efficiency measurement. It 
is additionally necessary to account for uncertainty in the amount of material in the inner detector region, which also 
affects the absolute track-finding efficiency. We make a conservative estimate of the systematic uncertainty by varying 
the total amount of material by 25%. This rather large variation in the amount of passive material in the detector 
translates into a systematic uncertainty on the tracking efficiency of approximately 2.6%. 
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J. Efficiency Summary 

For a given signal event that is within the acceptance of the detector, the absolute efficiency to reconstruct the 
event is the product of the efficiency terms outlined in this section. In many cases, the efficiency depends upon the 
track momentum as well as its position within the detector. The details of these efficiency parameterizations can be 
found in Ref. [2^ . For signal events that are well away from turn-on effects, the overall trigger and reconstruction 
efficiency is approximately 35%. 

VII. ACCEPTANCE 

Our detector acceptance is defined by the selection criteria and the active regions of the detector that cover those 
criteria. We determine the detector acceptance from MC simulation by ffi'st generating single b hadrons according to 
measured momentum and rapidity spectra. These hadrons are then decayed using the evtgen package [30] and then 
fed into a GEANT simulation of the detector. The MC simulation of the detector is utilized to map out the active 
regions of the detector. At this stage, we are not using the simulation to quantify detector inefficiencies, since those 
are measured using techniques described in the previous section. We generate b hadrons based upon two distinct 
distributions in transverse momentum and rapidity. The first is based on NLO [l|, 0| QCD theory with MRSDq 
parton distribution function (sij . a b quark mass of 4.75GeV/c^, and renormalization and factorization scales equal 
to ^ + py. Peterson fragmentation is used, with ep — 0.006. The other spectrum is based on the measured 
Hb J/ipX cross-section [3], which we refer to below as the "CDF MC sample." We use the production fractions 
fu-fd- fs of 1:1:0.270 [13]. We assume that Af, decays provide a negligible contribution to our signal because of the 
low probability for a semleptonic 6-baryon decay that would also yield a D'^ meson. 

To demonstrate the degree to which the MC simulation reproduces the data, Figs.[5l^a,b,c) show the pr distributions 
of the jjL" , /i, and for the MC simulation after all selection criteria and efficiencies have been applied, and for 
data after all selection criteria and corrections for the SVT efficiency. The data are seen to have higher prif^^D'-') than 
the MC simulation based upon the Hi, J/ipX cross-section result 8], while they are seen to have lower pt{^J.~ D^) 
than the spectrum from the MRSDq MC simulation, referred to below as the "MRSDq MC sample." We therefore 
measure the cross sections by taking the average of the cross sections from the two methods and assign one-half the 
difference as a systematic uncertainty. 

We perform additional comparisons between the data and the two MC simulated samples. Distributions such as 
r}{n~D°), pT{fJ.~), and pt{D^) are all seen to agree well between the data and both MC simulated samples. These 
comparisons can be found in Ref. [IB] • Table [ill shows the acceptance times efficiency (excluding the SVT efficiency) 
as a function of pxiHb) bin calculated using the CDF and MRSDq samples [26]. 

There is significant uncertainty in the branching ratios for iJ;, — > /i~ ^ D**X where the D** decays to D^X . An 
incorrect modelling of this fraction can affect the pT{pi^ D'^) spectrum and therefore the acceptance. To conservatively 
assess the systematic uncertainty on the measured cross section, we modify B{Hb fi^D**) by ±50%. 

VIII. BACKGROUNDS 

While fitting the D'^ peak to a Gaussian allows us to determine how many D'^ events are in our sample and account 
for noncharm background, it is possible that there are real events that are not part of our Hi, — > ^~ signal. 

One source of this background is direct cc production. It is possible for direct charm to mimic our signal by having 
one charm quark decay to a muon, and the other charm quark fragment into a . The second background source is 
bb production, where one b hadron decays to D'^X, and the other b hadron follows the decay chain H^ D ji^X. 
Since these background sources are irreducible, we estimate their contributions and correct the total number of 
observed signal events [1^. 

A. cc Background Estimate 

In order to estimate the fraction of events in our sample from direct charm, we use the impact parameter of the D'^ 
[28| . As the D'^ mesons from direct charm are created at the primary vertex, they should point back to the primary 
vertex within an error determined by the detector resolution. The from b decays, on the other hand, are created at 
the secondary b vertex, and are less likely to point back at the primary vertex. The D'^ impact parameter distribution 
for our signal sample is shown in Fig. [6l Although dominated by from &-hadron decays, this distribution contains 
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a prompt component arising from direct cc production. We use MC simulation to determine the expected impact 
parameter distributions for D'^ mesons from 6-hadron production and direct production. 

To determine the impact parameter distribution for D'^ mesons from b decays, we use the same MC simulated 
sample that is used to find the acceptance above. In order to get the do distribution of the from direct charm, 
we generate D'^ — > K~it'^ events originating from the primary vertex. The generated prompt D'^ K~'it'^ sample 
has lower pt{D^) than does our signal sample, so we reweight the px spectrum of the MC simulated so that it 
matches the pt{D^) spectrum in Hb fj,~D^ events from the pythia MC simulation [s^. 

In the direct charm and Hi, samples, we calculate the impact parameter using generator level MC simulated 

quantities. In order to estimate the uncertainty on the Z)° impact parameter due to detector resolution, we look at the 
sideband subtracted impact parameter uncertainty distribution in our signal sample acquired on the /xSVT trigger. 
We see that the average uncertainty is approximately 34 /xm, so we smear the MC simulated impact parameters using 
a Gaussian resolution function with a standard deviation of 34 /i m. We do a bin-by-bin fit to determine the ratio 
of events from direct charm to those from b hadron decays. We find that using the c?o distribution from the Hf, events 
generated with the pr spectrum measured in the CDF Hb J/iJjX cross section measurement we measure a 
charm fraction of (6.3 ± 2.1)%, while using the do distribution from the b events in the MRSDo sample gives a charm 
fraction of (5.4 ± 2.1)%. So as a final result, we use the average charm fraction of (5.9 ± 2.1(stat) ± 0.4(syst))%. 

For the differential cross section, we must also account for the pt spectrum of /i^D'^ candidates arising from direct 
cc production. For normalization, we use the total cc yield as calculated above. To determine the pt{D^) spectrum, 
we use MC simulation, where the generated events are reweighted to match the predicted direct charm spectrum [IJl . 
From this reweighted MC sample, we estimate the bin-by-bin contribution from direct cc background. 

B. bb Background Estimate 

The presence of bb events can mimic the signal in the cases where the /i^ and D'^ come from different b hadrons. An 
example of this is where the b hadron decays as 6 ^ D'^X and the other follows the decay chain b c ^ n~X. This 
configuration provides a ^~ candidate with the proper charge correlation, constituting an irreducible background 
to genuine H^ jj.^D'^ decays. 

To investigate this background source, we look at events where one b hadron decays to a Z?*^ meson, and the other 
follows the decay chain b fJ-^X , providing the wrong-sign ji^ D'^ combination. The right-sign hh background is then 
estimated by counting wrong-sign /i+Z?" events and correcting this by the expected ratio of the two decay chains 
b jX^ and 6 ^ c ^ . 

We begin by looking for wrong-sign jj.^D'^ events, where we have a K^n^ peak in events where the pion 

and muon have the same charge, = q^r- This is complicated by two factors. First, there is a large background 
in the wrong-sign sample from real signal events where the K~'!t~^ decay is incorrectly reconstructed 

as -D" K^n^ . The reflection does not yield a narrow D° peak in the pion-kaon invariant mass distribution but 
instead yields a broad peak which we parameterize using MC simulated decays. We use this parameterization and 
flx normalization of the reflection peak relative to the signal peak based upon the MC simulation. The MC simulated 
wrong-sign mass distribution is found to be insensitive to the input pt(&) spectrum. 

The second complication is that the doubly Cabibbo-suppressed decay of the K'^tt~ is also expected to 

yield a small peak in the K-n invariant mass distribution. Based upon measured branching ratios j3^ . the number of 
doubly Cabibbo-suppressed events is expected to be 0.0036 ± 0.0003 times the number of events in our signal peak. 
After correcting for these effects, we flnd a wrong-sign peak that is (4.4 ± 1.7)% the size of the right-sign peak. 

To convert this number into an estimate of the number of right-sign bb events, we look at a generator level sample 
of bb MC simulation generated using PYTHIA and decayed with evtgen using our best knowledge of branching 
ratios. In this sample, we look for events where one b decayed to a muon, and the other one produced a In 
order to get a sufficiently large event sample, we use somewhat looser selection criteria on the bb than are in our 
sample, and do not apply any of our efficiency measurements to the events. We flnd the ratio to be N(b — > c — > 
li~)/N(b fi~') — 0.23 ± 0.02 (stat). To take into account the effect of the looser event selection criteria, we apply 
a 50% systematic, giving a final ratio of 0.23 ± 0.02(stat) ± O.ll(syst). Combining this result with the wrong-sign 
contribution of (4.4 ± 1.7)% yields a bb background fraction of (1.0 ± 0.40 ± 0.48)%. 

C. Hb DD and Hb Dt Backgrounds 

In addition to backgrounds from cc and bb events, a background also comes from the decays of a single b hadron. 
These events can occur when the decay b —^ ccs or b —^ ccd is followed by one of the charm quarks decaying to a 
muon, and the other decaying through D° K^tr^ . These events also can come from a decay of 6 — > ct~T>t, where 
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the tau lepton decays to a muon, and the c decays through K~Tr~^. Because these events come from real Hi, 

decays, the impact parameter of the D° will not necessarily be consistent with production at the primary vertex. 
The muon and the pion from the decay will have opposite charges, as is the case in b ^ c^i^V^ events. We account 
for this contribution by retaining events with Hf, —> DD and iJf, — *■ Dt decays in the MC simulation utilized for 
the acceptance measurement and correct for it in the unfolding procedure described in Sect. HXl Since the branching 
fractions for some of these decays are not well known, we allow the number of events from Hi, — > DD and H}, Dt 
to vary by ±50%, which translates in a 2.5% systematic uncertainty on the measured cross section. 



D. Fake Muons 



In fitting the I?" mass distribution, we account for combinatoric background using the sidebands of the mass 
distribution. This technique is not available for hadrons that are misidentified as muons, so we must estimate the 
contribution of fake muons independently. 

For cc background events described above, the background is assessed by looking at the impact parameter of the 
-D", independent of whether the muon is real or fake. Therefore, the only source of fake muons that we must account 
for are fake muons reconstructed in association with real Hi, D^X decays. To evaluate this, we use a sample of 
£)+ K~TT^TT~^ decays from the heavy flavor trigger path that is enhanced in Hi, D~^X decays. We measure the 
rate of tracks near the _D+ which are fiducial to the CMU and CMP detectors and have pt > 3GeV/c. Because 
of their smaller nuclear interaction cross section, kaons are misidentified as muons at a higher rate than pions. To 
conservatively assess a muon misidcntification rate, we assume that all of these tracks are kaons. We then apply a 
0.5% probability that a fiducial kaon is reconstructed as a muon, as measured in Ref. [s^. From this, we determine 
that the rate of fake muons in our signal sample is approximately 0.05%, small enough to neglect. 



IX. Ht fj,-D''X CROSS SECTION 

In order to extract the 6-hadron differential cross section from the measured pt{h~ D'^) distribution, we unfold the 
observed pt{Hi,) distribution on the basis of the MC simulated signal. We determine a weight Wij that is defined as 

where N{[Hi,]^ — > [pL^D^]-) is the number of b hadrons with \y\ < 0.6 in pt{Hi,) bin i decaying to /x~Z3° in pt{^^ D^) 
bin J, which pass all selection criteria, and N{[^~ D'^]^) is the number of D^ events in pT{^ji~ D'^) bin j, which pass 
all selection criteria. 

We generate MC simulated Hi, events with \y{Hi,)\ < 0.8 as described in Section [VIII to properly model the 
detector acceptance. We also include events where the ^~D° is produced via a Hi, ^ DD or Hi, Dt decay in 
the denominator in order to take that background into account. For the determination of the weights Wij, it is the 
shape of the pt{Hi,) distribution, not its absolute normalization, that matters. We list the weights from the CDF 
and MRSDq MC simulated samples in Tables IIIII and IIVI Then, to unfold the number of b hadrons N^'' in a given 
PT{Hh) bin i we use the formula 



N 



(4) 



where Nj' ^ is the number of events in pt(m bin j from data, shown in Table IH We now have all of the terms 
from Eq. ([1]) required to extract the cross section. 

The statistical uncertainty on the cross section for pt{Hi,) bin i is given by 



Hb\2 



N 

E 



(5) 
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where and are the acceptance and efhciency for the bin i respectively. The statistical uncertainty on the total 
cross section is obtained by summing the uncertainties in quadrature over all bins of pt{Hi,). 

We check the analytical calculation of the statistical uncertainty using MC simulated pseudo experiments which 
model the signal and background distributions along with bin migration of the signal. We verify that the MC 
simulation yields the same statistical uncertainty on the differential and integrated cross section as we observe using 
the above analytical formula. The uncertainty is calculated by generating 1000 distinct MC simulated experiments 
corresponding to the total sample /i^D" yield we observe. The distribution of fitted yields in each px bin is observed 
to be Gaussian, and we take the standard deviation of the distribution to be the statistical uncertainty in each bin. 

The systematic uncertainties are determined in general by varying the efficiency or acceptance in question by 
±1(T and calculating the shift in the measured cross section. In the case of the efficiency measurements, which are 
parameterized based upon fits to the data {e.g. , the level 1 trigger efficiency is parameterized as a function of muon 
Pt and 77) the variation in efficiency is taken from the error matrix determined in the fits. 

Using this technique, the uncertainty derived based upon the SVT efficiency fit appears to be underestimated. The 
variation in the measured points is not reflected by the statistical uncertainty in the fit for the SVT efficiency. In this 
case, we assess the systematic uncertainty by replacing the fitted efficiency with the bin-by-bin points that are derived 
in the efficiency measurement. The difference between the result and the value obtained from the SVT efficiency 
data points is taken as a fit modelling uncertainty and combined in quadrature with the la systematic uncertainty 
described above to obtain the total systematic uncertainty ^26]. 

The systematic uncertainty on the MC simulated pt spectrum is estimated by taking the fractional difference 
between the cross section found using the Hi, J/ipX and the reweighted MRSDq NLO MC simulated sam ples . 
The degree to which final state radiation affects the lineshape is assessed with the PHOTOS MC simulation [3q . 
Systematic uncertainty contributions are shown in Table fVl It is worth noting that the systematic uncertainty arising 
from the MC simulated pT shape on the total cross section is less than the systematic uncertainty due to the MC 
simulation px shape in most of the px bins. This is because uncertainty on the pt shape primarily affects the fraction 
of the total cross section in a specific pr bin. When integrating over all bins, this effect is diluted. Also, when finding 
the acceptance for a specific px bin, one can only use the MC simulated events within that bin, while the overall 
measurement uses all of the MC simulated events. This leads to a lower statistical uncertainty in the MC simulated 
samples. 

After applying all corrections to the data, we obtain the total cross section times branching ratio 

a{pp HbX) X B{Hb ^I'D^X) x B{D° K-n+) ^ 3.46 ± 0.14(stat);o 38(syst) nb 

for b hadrons with Pt> 9GeV/c and \y\ < 0.6. Using the current world average branching ratios B{Hi, fjL^ D^) = 
(6.84 ± 0.35)% [Iltsl and B{D^ K-tt+) = (3.89 ± 0.05)% [H,!!!, we measure a total cross section of 

a{pp HbX) ^ 1.30 ± 0.05(stat) ± 0.14(syst) ± 0.07(6) ^lh 

for h hadrons with px > 9GcV/c and \y\ < 0.6. The differential cross section times branching ratio is shown in 
Table [VTl and displayed in Fig. [Tfa). Figure [Tl^b) shows the differential cross section. 

X. Hb fi-D*+X CROSS SECTION 

For the Hb fi^D*^, D*^ — > Z3'^7r+, D*^ K^ir^ channel, the efficiencies are the same as those for the 
Hb ^JL^ , D'^ if ^7r+, except that we must additionally account for the soft pion identification efficiency. 

The systematics for the D*^ measurement are the same as those for the D'^ measurement listed in Table fVl unless 
mentioned specifically. 

After applying all corrections to the data, we obtain the total cross section times branching ratio Hb —f fi^D*^ , 
with D*+ D°'K+ and D° K~'k+ decays of 

C7{pp HbX) X B{Hb ^ fi-K-n+TT+X) = 1.05 ± 0.08(stat)+°;|5(syst) nb, 

for b hadrons with pT > 9GeV/c and \y\ < 0.6 and using the notation B{Hb fi^ K^ti^tt^ X) as shorthand 
for the product of B{Hb n^D*^), B{D*~^ D^tt^) and B{D^ K^ir^). Using the current world average 
branching ratios B{Hb -> M"^*^I.^J|-75 ± 0.19)% |li, HI, B{D*+ D^tt+) = (67.7 ± 0.5)% [H HI, and 
B{D° K-n+) = (3.89 ± 0.05)% [H Hi , we find a total cross section of 
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(t{pp ^ HbX) = 1.45 ± 0.11(stat)Io;2i(syst) ± 0.10(B) 

for b hadrons with pT > 9GeV/c and |y| < 0.6. The differential cross section times branching ratio is shown in 
Table IVIIl and displayed in Fig. [5] We also show the differential cross section using branching ratios measured 
elsewhere. 

We consider the Ht j^r D*^ X cross-section results to be a modest extension and cross check of the primary 
Hjj ^~ D^X cross section result presented in this paper. The ^~ D*^ X sample is of limited statistical power and, 
being a subset of the ii~ D^X sample, there is nothing to be gained by averaging the two results. 

XI. COMPARISON WITH THEORY AND PREVIOUS MEASUREMENTS 

Our results agree well with previous measurements of the Hi, cross section at CDF [s^ and with the predictions of 
FONLL [l3|. The CDF collaboration previously measured the cross section using — > J /ifjX of 1.34 ± 0.07/ib 
in the same range of pt and y [1]. The FONLL prediction for the total cross section with pxib) > 9.0GeV/c and 
\y\ <0.6 is 1.39t[J:gMb 0. Figure [9] shows the differential cross sections plotted together. 

To compare these results with the cross section measured in the exclusive decay channel 5+ J/ipK'^ [l3|, we 
must account for the different rapidity ranges. Reference Jl 2] measured the cross section for \y\ < 1, while the result 
presented here and the inclusive J/ipX result from Ref. [3| measure the cross section for \y\ < 0.6. Figure [TOl shows 
the comparison between these results, where we have scaled all cross sections to |y| < 1 assuming that inclusive Hh 
production is flat in rapidity. Furthermore, we correct for the differences in fragmentation fractions for (/„) and 
Hb ifu + fd + fs) [13 • Figure (TU] shows good agreement between these recent 5-hadron cross section measurements 
at -^s — 1.96 TeV, as well as agreement between these complementary measurements and the FONLL prediction. 

The results presented here also compare favorably to the recent theoretical work using the general-mass variable- 
flavor- number scheme (GM-VFNS) presented in Ref. fssj. 

The trigger and analysis strategy shown here for the /i~Z?° and p^D*^ final state are quite different than those 
utilized in measurements using J ftp final states. In addition to being statistically independent, many of the systematic 
uncertainties of this result are distinct from those of previous measurements. This result therefore provides a unique 
and independent measurement of the Hf, cross section at the Tevatron. 

XII. CONCLUSION 

We have measured the 5-hadron production cross section in iJf, —> fi^D'^X and Hb ^jT D*^ X decays using 
the CDF II detector at the Fermilab Tevatron. The muon plus charm hadron decay signature was acquired using a 
dedicated trigger that takes advantage of the muon signature as well as the long Hb lifetime. By selecting an event 
sample based upon decay length information measured in the trigger, it is possible to acquire clean Hb signal samples 
at lower px^Hb) than have been measured before. Complementary data samples are used to measure the trigger 
efficiencies for this analysis and Monte Carlo simulation is used to measure the acceptance. After accounting for 
backgrounds that mostly arise from heavy flavor [bb and cc production) the cross section is measured to be 

(t{pP HbX) = 1.30 ± 0.05(stat) ± 0.14(syst) ± 0.07{B) ph 

for b hadrons with pT > 9 GeV/c and \y\ < 0.6. We additionally present the differential cross section da/dpr in this 
kinematic region, which extends to lower pxiHb) than previous measurements in semileptonic modes. We find the 
results presented here to be in good agreement with other recent measurements of the Hb cross section and in good 
agreement with fixed order next-to-leading logarithm calculations. 
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TABLE I: The signal yields for the number of D" and D*~^ events per pt{ii~ D") bin before (and after) correcting for SVT 
efficiency. Since the SVT efficiency depends upon the relative track isolation, a correction is performed based upon the topology 
of the signal events. 

pril^'D"^) [GeV/c] D'J signal yield D* signal yield 

Measured Corrected Measured Corrected 

9Ai 867.9 ± 53.5 1040.6 ± 64.1 82.3 ± 10.5 96.1 ± 12.2 
11-13 863.1 ± 45.8 1034.4 ± 54.9 142.6 ± 13.2 170.0 ± 15.7 
13-17 1016.8 ± 46.3 1192.4 ± 54.3 236.4 ± 16.9 283.4 ± 20.2 
17-29 669.6 ± 38.3 781.9 ± 44.7 169.6 ± 14.1 199.8 ± 16.6 
29-40 67.0 ± 12.2 80.1 ± 14.6 14.3 ± 4.2 18.7 ± 5.6 



TABLE II: The acceptances times efficiencies using the CDF phenomenological spectrum (CDF Sample) and MRSDo NLO 
(MRSDo Sample) MC simulation. The average priHt) for the two spectra vary by no more than 1% in any transverse 
momentum bin. 

PT{Hb) [GeV/c] Average pT(-ff6) [GeV/c] Overall efBciency 

CDF Sample MRSDp Sample 

9^Tl 9^9 1.47 X IQ-'^ 1.51 x 10"^ 

11-13 11.9 6.12 X 10-3 6.16 x lO'^ 

13-17 14.8 1.38 X 10-2 1.40 x 10"^ 

17-23 19.6 2.68 x lO'^ 2.70 x lO'^ 

23-29 25.5 3.90 x lO'^ 3.84 x 10"^ 

29-40 33.2 4.61 x 10"^ 4.63 x 10"^ 



TABLE III: The weights Wij from the CDF phenomenological spectrum MC simulation with all momenta in GeV/c. 
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TABLE IV: The weights Wij from the NDE+MRSDo MC simulation with all momenta in GeV/c. 
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-4 


2.28 


X 


10- 


-3 


1.27 


X 


10- 


-2 


1.91 


X 


10- 


-1 


6.00 


X 


10- 


-1 


PT{Hb) > 40GeV/c 


< 1 


X 


10- 


-6 


7.90 


X 


10- 


-5 


7.30 


X 


10- 


-4 


1.76 


X 


10- 


-2 


3.67 


X 


10- 


-1 



TABLE V: The systematic uncertainties of the da{pp — » Hb)/dpT measurement. All numbers listed in percent. 



Pt range [GeV/c] 



Source 


All 


9-11 


11-13 


13-17 


17-23 


23-29 


29-40 


Luminosity 


6.0 


6.0 


6.0 


6.0 


6.0 


6.0 


6.0 


LI Efficiency 


2.6 


2.6 


2.6 


2.6 


2.5 


2.5 


2.5 


L2 XFT Efficiency 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


L2 SVT Efficiency 


1.3 


1.6 


1.2 


0.8 


0.7 


0.7 


0.9 


L3 Efficiency 


0.2 


0.2 


0.2 


0.2 


0.2 


0.1 


0.1 


CMU Efficiency 


1.6 


1.6 


1.6 


1.6 


1.6 


1.6 


1.6 


CMP Efficiency 


-1-1.4 
-1.3 


+ 1.4 
-1.3 


+ 1.4 
-1.3 


+ 1.4 
-1.3 


1.4 


1.5 


1.7 


SVX II Efficiency 


-1-1.5 
-1.0 


+ 1.5 
-1.0 


+1.5 
-1.0 


+1.5 
-1.0 


+1.5 
-1.1 


+ 1.5 
-1.1 


+ 1.5 
-1.0 


SVX II 2nd track Efficiency 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


COTEfF. (DO) 


+.'i.i 

-4.1 


+3.2 
-4.3 


+3.3 
-4.2 


+3.0 
-3.8 


+2.9 
-3.5 


+2.8 
-3.3 


+2.7 
-3.1 


COT EfF. {ii-D*+) 


+ 7.Q 
-8.9 


+ 7.9 
-9.7 


+ 7.1 
-9.6 


+6.1 
-8.2 


+ 5.3 
-6.1 


+4.9 
-5.0 


+4.5 
-4.3 


Vertex + Fit Eff 


0.4 


0.4 


0.4 


0.4 


0.4 


0.4 


0.4 


cc Background 


4.7 


5.3 


5.8 


3.8 


2.5 


2.1 


1.7 


hb Background 


0.6 


0.6 


0.6 


0.6 


0.6 


0.6 


0.6 


rfo Smearing 


-1-0.9 
-2.6 


+1.1 

-2.6 


+0.8 
-2.8 


+0.7 
-2.6 


+0.8 
-2.3 


+ 1.0 
-2.5 


+ 1.5 

-1.4 


final state radiation 


1.2 


1.0 


1.8 


1.8 


0.7 


2.2 


4.0 


B{Hi ii-D**X) in-D") 


1.9 


2.1 


1.7 


1.6 


2.0 


1.5 


1.7 


B{Hb ii-D**X) {ii-D*+) 


1.9 


2.0 


1.9 


1.8 


2.3 


1.7 


2.0 


B{Hb -> DD), B{Hb ^ Dt) 


2.4 


2.4 


2.4 


2.5 


2.6 


2.3 


2.4 


MC Stat. ifJ.-D°) 


0.5 


1.7 


1.1 


0.9 


1.0 


1.6 


2.1 


MC Stat. {ii-D* + ) 


1.2 


5.8 


3.4 


2.2 


2.2 


3.3 


4.4 


MC Pt shape {n'D^) 


3.4 


5.7 


3.0 


1.7 


2.2 


2.0 


2.5 


MC Pt sliape {i^-D*+) 


4.4 


5.7 


4.4 


3.9 


4.0 


0.6 


2.8 


Total {fi-D°) 


+ 10.4 
-11.0 


+ 11.9 
-12.4 


+ 11.0 

-11.5 


+9.6 
-10.1 


+9.2 
-9.6 


+9.2 
-9.6 


+ 10.0 
-10.1 


Total {fi-D*+) 


+ 12.5 
-13.9 


+ 14.9 
-16.0 


+ 13.5 
-15.1 


+ 11.7 
-13.1 


+ 11.0 
-11.5 


+ 10.3 
-10.5 


+ 11.5 
-11.3 
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TABLE VI: Differential cross section dcr{pp — > Hf,)/dpT and differential cross section times branching fraction of Hi, — » D°X, 
K~n'^ , where the first uncertainty is statistical, and the second is systematic. In the case of the differential cross section, 
the third uncertainty arises from the applied branching fractions. 



PT bin [GeV/c] da/dpr x B [pb/(GeV/c)] da/dpr [pb/(GeV/c)] 

9^11 762 ± 50tll 287Tl9t|f±15 

11-13 403±18t^^ 152±7±17±8 

13-17 179±6ti^ 67.3 ± 2.3^J?;8 ± 3.6 

17-23 49.7 ±1.5^^1 17.8 ± 0.6tJ] ± 0.9 

23-29 13.1 ± 0.6^1 j 4.93 ±0.231^:4^ ±0.26 

29-40 3.43 ± 0.20tg j5 1.29 ± 0.08 ± 0.13 ± 0.07 



TABLE VII: Differential cross section da{pp — » Hb)/dpT and differential cross section times branching fraction of Hb — » 
fx~ D*'^X, £)*"*" — » D^tt"*", D° — » if "tt"*", where the first uncertainty is statistical, and the second is systematic. In the case of 
the differential cross section, the third uncertainty arises from the applied branching fractions. 

PT bin [GeV/c] da/dpr x B [pb/(GeV/c)] da/dpr [pb/(GeV/c)] 



9-11 226 ±301^6 312 ± 41+1;, ± 22 

11-13 122±10tJ| 168 ± 141:^^ ± 12 

13-17 56.0 ± 3.0t^ j 77.3 ± A.ltlo,^ ± 5.5 

17-23 15A±0.7tli 21.3 ±1.0l2:i± 1-5 

23-29 3.84±0.25t°;^° 5.30 ± 0.35 ± 0.55 ± 0.37 

29-40 0.95±0.09t^;iJ 1.31 ± 0.12 ± 0.15 ± 0.09 
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FIG. 1: The tv~^K~ mass for difFerent ranges of pt{h~ K~n^). Points indicate events where the muon and pion have opposite 
charges, which is the right-sign correlation, and the dashed histograms contain events where the muon and pion have the same 
charge, which is the wrong-sign correlation. The solid line is a fit to the right-sign data using a Gaussian plus linear background. 
The per degree of freedom for these fits range from 0.9 to 1.2. Because of low statistics, the yield in the highest pt bin is 
measured by sideband subtraction. 
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Kti mass (GeV/c ) 

FIG. 2: The -k^ K~ mass for pT{fJ-~ K~ tt'^) > 9.0 GeV/c. The sohd histogram shows events with the right-sign correlation 
between muon and pion, and the dashed histogram shows events where the muon and pion have the wrong-sign correlation. 
The solid line is a fit to the right-sign data using a Gaussian plus linear background. 
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FIG. 3: The mass difference between TisojtT^'^ and iv'^ K~ for different ranges oi pr {tJ-~ K~ iv^) . All events shown are selected 
such that the muon and the pion from the candidate have opposite charge. Points show events where the muon and soft 
pion have opposite charge, which is the right-sign correlation, and the dashed histograms represent events where the muon and 
soft pion have the same charge, which is the wrong-sign correlation. The per degree of freedom for these fits ranges from 
1.1 to 1.3. Because of low statistics, the yield in the highest pr bin is measured by sideband subtraction. 
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FIG. 4: Distribution of the difference between the ''^toft''^^^~ ^^i^ the -k^ K~ invariant masses for pT{ii~ K~ -k^) > 9.0 GeV/c. 
The solid histogram indicate events with the right-sign correlation between muon and pion, and the dashed histogram contain 
events where the muon and pion have the wrong-sign correlation. 
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FIG. 5: A comparison of the two MC simulated samples to data, plotted as a function of the pT of (a) the /i _D°, (b) the /i, 
and (c) the _D°. The data is seen to have a higher pt{^~ D'^) than the CDF MC sample, while the data is seen to have a lower 
PT{fJ-'D°) than the MRSDq NLO MC simulation. 
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FIG. 6: The D impact parameter distribution for the data after sideband subtraction. 
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FIG. 7: The differential cross section times branching ratio for Hb — + fj.~ D'^X, K~-k'^ is shown on the left, where BR is 

shorthand notation for the product of branching ratios, BR = B{Hi, /i^ D^X) x B{D° — > K~'k^). The uncertainties shown 
on each point include statistical and systematic uncertainties combined in quadrature. Incorporating measured branching ratios 
[3^ . the differential cross section is shown on the right. The uncertainties on each point include statistical, systematic and 
branching ratio uncertainties added in quadrature. 
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FIG. 8: The diflerential cross section times branching ratio for H^, fj, D*^X, D*"'" D'^'^toff ~^ ^ shown on the 

left, where _B_R is shorthand notation for the product of branching ratios, BR = B{Hb —* fi~ D'^X)xB{D'^ n'^)xB{D*'^ 
D^tt"''). The uncertainties shown on each point include statistical and systematic uncertainties combined in quadrature. Incor- 
porating measured branching ratios [s^l, the differential cross section is shown on the right. The uncertainties on each point 
include statistical, systematic and branching ratio uncertainties added in quadrature. 
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FIG. 9; The 6-hadron differential cross section for \y{Hb)\ < 0.6 from FONLL theory [lol ] compared with measurements from 
Ht n-D°X {D° K--K+) and Hi, fi-D'+X {D*+ D^-^toft and D° ^ K-tt+). 
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FIG. 10: Comparison of difFerential b cross-section results. The results presented are Ht J/ipX at ^/s — 1.96 TeV [1), 
B+ J/ipK+ at ^ = 1.96 TeV [l2|, and this result for Hi ^"D^X at = 1.96 TeV. All results are scaled to \rj\ < 1. 
For comparison, the FONLL [l^ calculation is included. 



